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Abstract 

We propose and analyze a variation of the Euler scheme for state 
constrained ordinary differential inclusions under weak assumptions 
on the right-hand side and the state constraints. Convergence results 
are given for the space-continuous and the space-discrete versions of 
this scheme, and a numerical example illustrates in which sense these 
limits have to be interpreted. 

Keywords: differential inclusion, state constraints, numerical method 
AMS subject classifications: 34A60, 65L20, 49J15 


1 Introduction 

Numerical methods for ordinary differential inclusions have mainly been con¬ 
sidered in the unconstrained case and for convex-valued right-hand sides. 
Summaries of early results can be found in the surveys [10] and [T3| . An ap¬ 
proach using subdivision and continuation techniques for time-independent 
systems with affine controls was given in [T8] . 

The variation of the Euler scheme proposed in [12] only uses extremal 
points of the right-hand side, so that a fully discretized scheme is obtained 
for polytope-valued right-hand sides. The high complexity of this method 
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motivated the detailed analysis of spatial discretization effects on the En- 
ler scheme given in [7]. A further reduction of computational costs can be 
achieved by tracking the boundaries of the reachable sets of the fully dis¬ 
cretized Euler scheme using only lower-dimensional data of the right-hand 
side, see [T6] . 

A very different type of algorithm has been proposed in [1], where a 
reachable set is computed by an optimal control routine that minimizes the 
distance to grid points in the relevant region at the terminal time. As the 
optimizer may hnd a local instead of a global minimum, the propriety of this 
method cannot be guaranteed. 

An implicit Euler scheme has been analyzed in [S], which is considerably 
more efficient than the explicit Euler scheme when applied to stiff differential 
inclusions. The disadvantage of this method is that it has to solve an alge¬ 
braic inclusion in every time step. This problem is successfully avoided by the 
semi-implicit Euler schemes discussed in [15], leading to a better performance 
of the algorithm. 

The Euler scheme for differential inclusions with nonconvex right-hand 
sides has been investigated in mi. Error estimates for the Euler scheme 
applied to a differential inclusion with convex, but one-sided Lipschitz right- 
hand side have been published in [9]. 

Up to our knowledge, the only convergence analysis of the Euler scheme 
for ordinary differential inclusions with state constraints has been given in 
the paper [5|. The analysis is based on stability theorems that quantify the 
relationship between the solution set of the constrained problem and the same 
differential inclusion without constraints. These results play an important 
role in optimal control theory and are still a subject of intensive research, see 
0.0 and the references therein. 

The drawback of these stability theorems in the present context are the 
strict assumptions on the compatibility between right-hand side and con¬ 
straints that are needed to obtain such a result. In applications like collision 
avoidance, see nn, na and the references therein, such assumptions do cer¬ 
tainly not hold. 

The present article is motivated by the lack of a convergence analysis for 
the Euler scheme in the present of state constraints under weaker conditions. 
In Section [2] we £x the notation, give a formal problem statement and for¬ 
mulate our basic assumptions, that the initial set is compact, the right-hand 
side is Lipschitz, the state constraints are upper semicontinuous in time and 
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that all data are compatible in the very weak sense that there is a nontrivial 
time interval on which a solution of the state constrained inclusion exists. 

Under these conditions, we prove some auxiliary results in Section |3l 
which are used in Section 0] to show that the solution sets of the inflated 
Euler scheme converge towards the exact solution sets in Hausdorff distance. 
In Section 0 we show that the same holds for a spatially discretized version of 
the scheme provided all constraints are relaxed in a proper way. A numerical 
example given in Section [ 6 ] illustrates in which sense this limit has to be 
interpreted. Note that convergence results for the inflated Euler scheme 
without spatial discretization are of interest for first-discretize-then-optimize 
strategies in optimal control and algorithms such as the one proposed in 
[1], which consider the Euler scheme as a constraint of a finite-dimensional 
optimization problem. 

At the end of the paper, in Section [71 we give a brief explanation how sta¬ 
bility results in the spirit of [5] and | 6 ] automatically yield linear convergence 
of the inflated Euler scheme in Hausdorff distance. 


2 Setting and notation 

Throughout this paper, the Euclidean norm and the modulus will be denoted 
by II ■ II : ^ R 4 . and | ■ | : R —R+. As some of the following concepts will 

be used not only in the state space R'^, but also in the space of continuous 
functions, we introduce them in a general Banach space E. 

Let A,R C E be compact. Then 

dist(A, B) := max min ||a — 6 || 

a&A b£B ' ' 

is called the Hausdorff semidistance between A and B. In the particular case, 
where A = {a} with a ^ E, the Hausdorff semidistance is the usual distance 
dist(a. A) from a point to a set, and we denote 

Bji{A) ■= {x E E dist(x. A) < R}. 

The Hausdorff distance 

dist_H-(A, B) := max{dist(A, B), dist(i?. A)} 

is a metric on the compact subsets of E. The metric projection from x ^ E 
to A is the set 

proj(a:. A) := {a E A : ||a; — a|| = dist(a;. A)}. 
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Since A is compact, the metric projection is nonempty. If, in addition, B is 
a Hilbert space and A is convex, then proj (x, A) is a singleton. 

By hh^’^([0, r], R'^) with r G (0, cxo), we denote the Banach space of abso¬ 
lutely continuous functions from [ 0 , r] to R'^, i.e. the space of all functions, 
which possess a weak derivative in h^((0, r), R'^). 

Let F : [0,T] x R'^ =1 R'^ be a set-valued mapping. A solution to the 
ordinary differential inclusion 


x(t) e F(t,x(t)) ( 1 ) 

is an absolutely continuous function x(-) G kh^’^([0, T], R'^) satisfying ([T]) 
almost everywhere in [0,T]. Depending on the application, solutions may be 
required to satisfy an initial condition 

x(0) G Xq (2) 

with initial set Xq C R'^ and state constraints 

x{t) G A{t) (3) 

with A : [0,T] =1 R"^. For r G (0,T], we denote the solution sets of the 
unconstrained and the constrained initial value problems by 

S'“(r) :={a:(-) G kF^’^([0, r], R"^) ; ([1]) holds a.e. in (0,r), (|2]) holds}, 
:={x(-) G 5'“(r) : x{t) satishes ([3]) for all t G [0,r]}. 

These sets will be approximated by the corresponding solution sets of the 
inflated Euler scheme. For G N \ {0}, set h]\f := T/N, dehne a temporal 
grid t 7 v,n = nh^, n = 0,... ,N, and let (3^, 6^ G [0, cxo) be numbers to be 
specihed later. Then any sequence {yN,n)n C R'^ satisfying 

yN,n+l £ yN,n + hi^F{ti^^n, yN,n) + (4) 

is called a solution of the inflated Euler scheme with step-size h^, which may 
as well be subject to a possibly relaxed initial condition 

dist{yN,o,Xo) < 5 n (5) 

and possibly relaxed state constraints 

dist(|/Ar,„, A(t7v,n)) < Sn- (6) 
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For u E {0,..., N}, we denote the solution sets of the discrete unconstrained 
and the discrete constrained problems with time horizon u by 

Sn{^,/3n,Sn) ■={{yN,n)n=o C : (SI) holds for 0 < n < - 1, dS]) holds}, 

S%{u,I3n,Sn) ■.={{yN,n)n=o ^ , ^n) ■ (E]) holds for 0 < n < p}, 

and we identify these sequences with continuous functions |/ 7 v(-) : [0,ph] -E 
E'^ via piecewise linear interpolation of the data {tN,n,yN,n)n=oy where the 
representant of the derivative is chosen piecewise constant on intervals of the 
form ^Af,n+l)• 

Throughout this paper, we posit the following assumptions on the ordi¬ 
nary differential inclusion and the constraints. 

Hypothesis A. The right-hand side F, the constraints A and the initial set 
Xq satisfy the following conditions: 

(1) The mapping F : [0, T] x E'^ ^ E"^ has convex and compact values and 
is L-Lipschitz, i.e. there exists some L > 0 such that 

dist(F(t, x), F{t', x')) < L\t—t'\-\-L\\x—x'\\ for all t, t' G [0, T], x, x' G E'^. 

(2) The state constraints A : [0,T] E'^ have closed values and are upper 

semicontinuous in the sense that t = lim„^.oot„, Xn G A{tn) and x = 
Wmn^ooXn imply x G A{t). 

(3) The set Xq C E'^ is compact, and the data F, A and Xq are compatible 
in the sense that there exists r G (0,T] with ^ 0. 

We will also posit assumptions on the blowup size /Sv of the Euler scheme 
and the relaxation parameter for the initial condition and the constraints. 

Hypothesis B. The parameters (3^ and 6 n satisfy the following conditions: 

(1) The blowup size satisfies 0 as N ^ oo. 

(2) The relaxation parameter satisfies 6 n ^ 0 as N ^ oo. 

The relaxation parameter dpf will not be needed before Section O We 
introduce it at this early stage, because we wish to provide a-priori bounds 
in Lemmadl which apply to all exact and numerical trajectories in this paper. 
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3 Auxiliary results 


Since S^{t) C *S'“(r) and (3n,Sn) C S'^{h', (3 n,^n), the bounds given 

below hold for the constrained problems as well. 

Lemma 1. We have the a-priori bounds 

dist(a:(t), Xo) < R{t) := — l)(sup ||F(0,a;)|| + 1) Vt G [0, r] 

xGXq 

for all x{-) e S'“(r) and 

dist{yN,n,Xo) < RNin) := sup ||F(0, x)||+l) Vn G {0,..., i/} 

xeXo 

for all {yN,n)n=o ^ f^N,dN)- Moreover, we have 

R ;= max{l, sup ||F(0, a:)|| + sup {dx + -^f-) + 1) < oo, 

xeXo vgn ^ 

we have R{t) < R for all t G [0, T] and Rxin) < R for all n ^ {0,..., N} 
and iV G N, and for all (t,x) G [0,T] x Bji{Xq) we have 

||F(t,x)|| <P:=L{T + R)+ sup ||F(0,a;')|| < oo. 

x'eXo 

Proof. The bounds for S^{t) and (3x, dx) follow by applying Gronwall’s 
Lemma and induction, respectively, in a straight-forward manner. Bounded¬ 
ness of F is a consequence of hypothesis (Al). □ 

The following lemma quantifies the error between exact trajectories and 
Euler trajectories with perturbed initial value. It will be used to determine 
an appropriate size (dx of the blowup for the inflated Euler scheme. 

Lemma 2. Let t G [0, T] and h > 0 be such that [t,t-\-h] C [0, T], let z G 
and let x{-) G W^'^{[t,t h], Br{Xo)) be a solution of ([T]). Then there exists 

V G F(t, z) such that 

\\x{t h) — {z -\- hv)\\ < (1 -|- Lh)\\x{t) — z\\ L(P l)h‘^. 


Proof. As the metric projection to a compact and convex set is continuous 
according to [2], Theorem 9.3.4], we may define 




V := 


pToi{x{s),F{t,z))ds, 
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and V G F{t, z) holds by [ISl Theorem 1.6.13]. Because of Lemma [H we have 

pt-\-h 


Ixfs') - x(t)\\ < 


for s G [t,t + h], and we hnd 






x{s)ds — hv\\ < 


\x{e)\\de < Ph 


|i:(s) - pToi{x{s),F{t,z))\\ds 




dist{x(s), F{t, z))ds < / dist{F{s,x{s)), F{t, z))ds 




< 


L{s — t) + L\\x{s) — x{t)\\ + L\\x{t) — z\\dt 


< Lh\\x{t) — z\\ + L{P + l)h^, 


so that 


ct-\-h 


\x{t P h) — [z + hv)\\ < llx(f) — .^11 + II y x{s)ds — hv\\ 

< (1 + Lh)\\x{t) - z|| + L{P + l)h^ 


□ 


The following result, which is proved by compactness type arguments, is 
the core of this paper. In the current situation, it replaces Filippov’s theorem 
(see [H Theorem 2.4.1] for a continuous and [21 Theorem 2.2] for a discrete 
version of this result for Lipschitz right-hand side), which yields estimates 
for the semi-distances between numerical end exact solutions in the absence 
of state constraints. Here and later on, we equip the space R x R'’* x with 
the norm | ■ | + || • || + || • ||- 

Lemma 3. Let (rAr)Ar C (0,T] be a sequence with r := liniAr^oo'rv > 0, let 
c(-) G L^(0,T) and let xn{-) G hF^’^([0, tat], R'^) he a sequence of functions 
satisfying 

||iv(^)|| < c(t) for almost every t G (0,rAr), (7) 

sup dist((f, XAr(t), i;Ar(f)), graph(F)) —)■ 0 as N ^ oo, ( 8 ) 

tG[0,Tjv] 

dist(a:Ar(0), Xo) — t 0 as N ^ oo, (9) 

sup dist{xN{tN,n), A{tN,n)) 0 as N ^ oo. (10) 
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Then there exists some x{-) G 5'‘^(r) and functions XAr(-) G r], R'^) 

that coincide with xn{-) on [ 0 , tat] such that 

sup ||xAf(^) — 2 ;(t)|| —0 as N ^ oo, 

4G[0,t] 

xn{-) i(-) in L^{{0, t), K,"^) 

along a subsequence. 

Proof. Define 

;;; t G [0, min{r 7 v, t}], 

^ ' \ xnItn), t G (min{rAr,r},r]. 

Then the sequence {xn{-))n C y4C([0, r], R'^) satisfies condition ([7j) with r 
instead of tn, and it inherits conditions (IHl), (Ej) and ffTOl) with min{rAr,r} 
instead of tn- Because of (171), the functions {xn{-))n are uniformly bounded, 
and [U Theorem 0.3.4] guarantees that there exists some x{-) G hT^’^([0, r], R'^) 
such that 


sup ||xAr(t) - a:(t)|| -)■ 0, (11) 

tG[0,T] 

xn{-) ^ x (-) in L^((0,r),R'^) 

along a subsequence. By (IH]), the convergence theorem [H Theorem 1.4.1] 
applies, and together with (E]) and compactness of Xq, we obtain x{-) G -S'“(r). 

Let us check that x{-) G For any t G [0, r), we have tArqt//ijvj < ^ < 

min{rAr,r} for almost all N and t = limjv^oo• By ffTOj) . the points 
zn ■= p^oj{xN{tN^it/hN\), A{tN,it/hN\) satisfy 

||5v(iv,p/hjvj) ~ ^v|| 

= dist{xN{tN,[t/hj^\), A(tw,Lt//ijvj)) 0 as X -)> CX). (12) 

By continuity of x{-) and because of statements ffTTl) and (El), we have 
||a:(t) - Zn\\ 

< ||a;(t) - x{tjs[^it/hM\)\\ + \\^{tN,lt/hNi) - ^v(^Af,Lt//xivj)ll + ll^w(^Ar,Lt//iivj) “ 

—)■ 0 as X —)■ oo, 

so that hypothesis (A2) implies 

x{t) = lim Zn e A(t). 

N^oo 
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Since t G [0, r) was arbitrary, we also have x{t) G A.{t) by continuity of x{-) 
and hypothesis (A2). Consequently, we have x{-) G S'^{t). □ 

As a consequence, we obtain the following statement about the life span 
of the solutions to the the state constrained differential inclusion. 

Corollary 4. The number 

T* := sup{r G (0,T] : S"'(r) ^ 0} 

is, in fact, a maximum. 

Proof. By hypothesis (A3), there exists r G (0,T] with S^{t) ^ 0, so that t* 
is well-dehned. By dehnition of the supremum, there exist X 7 v(-) G 
N E N, such that -E- r*. According to Lemma |3l they possess a limit, 
which is a member of □ 

4 The inflated Euler scheme 

The following example shows that it is necessary to modify the Euler scheme 
to treat state constrained problems. 

Example 5. Consider the constrained differential equation 

Xi =-X 2 , X 2 = XI, xl + xl = l, Xi < 0, a;i(0) = 0, X2(0) = 1, 

with unique solution 

Xi{t) = sin(t + 7r/2), X 2 {t) = cos{t + 'x jT), 

and maximal interval of existence [0,7r]. For any step-size h > 0, the Euler 
scheme maps (0,1) to the point (—h, 1), which does not satisfy the constraints, 
so that no feasible approximate trajectory exists. 

Motivated by Lemma [2] and the above example, we choose a blowup size 

I3N-.= L{P + I)hl. (13) 

In this context, there is no need to relax the initial condition or the state 
constraints, so that we have = 0. A relaxation will become necessary 
in Section |5l where not only the time interval, but also the phase space is 
discretized. 

The following proposition states that with this choice of /Sat, the exact 
dynamics are overapproximated by the inflated Euler scheme. In particular, 
Euler trajectories exist at least as long as exact trajectories. 
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Proposition 6. For any x{-) G 5'‘^(r) with r G (0,r*], there exists a sequence 
{yN)N withy^ G S'^([r//i7vJ,/3Ar,0) such that 

x{tN,n) = yN,n, Ti = 0, . . . , [t/ZiatJ . (14) 

Proof. The existence of a sequence {yN)N with fflTl) can be obtained by suc¬ 
cessively applying Lemma [2] and taking flT^ into account. □ 

The statement of Proposition [ 6 ] can be reformulated in terms of the Haus- 
dorff semidistance as 

sup inf max \\yN,n - x(tN,n)\\ = 0. 

x(-)eS^(T) yM^‘S^([T/hjvj,/3N,0) n=0,...,[T/h/^] 

The following proposition guarantees that the opposite Hausdorff semidis¬ 
tance between numerical and exact trajectories converges to zero. In partic¬ 
ular, Euler trajectories do not live longer asymptotically than exact trajec¬ 
tories. 

Proposition 7. For any r G (0,r*], we have 

sup inf max \\yN,n - x(tN,n)\\ ^ 0 (15) 

yjvGS^([T/fcjvJ,/3jv,0) 3;(-)eSqT) n=0,...,[T//ijvJ 

as N ^ oo. For any r G (r*,T], there exists E ¥l such that 


*^w(L^/M,/9iv,o) = 0 yN>Nr. 

Proof. Let r G (0, r*]. By Corollary 0] and because of Proposition [ 6 l we have 

and ^^(Lr/M,/5^,O)^0. 

Assume that statement flTHl) is false. Then there exist e > 0 and a sequence 
yN G S'^( [r//iArJ,/^AT, 0), iV G N' C N, such that 


inf max \\yN,n - x{tN,r. 


> e. 


(16) 


The linearly interpolated Euler trajectories |/Ar(-) satisfy the assumptions of 
Lemma [3l because we have tN,[T/hNi estimate 

||^v(s)||<P Vs G [0,tjvq,/A^j] 
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holds according to Lemma [H we have 

dist((s, 2 /iv(s), y 7 v(s)), graph(F)) 

^ IK®) 2 /v(s), yv(s)) — (iAr,[s//ijvJ) 2 /Af(^V,[s/h]vj))yw(iv,[s/hivj))ll 
+ dist((tAr,Ls//i^j, yN{tf^^is/hN\)^yN{tN,ls/hN\)yS^^VH^)) 

^ (^1. -\- P')h]\f -\- I3]\f / h}^ —^ 0 as N —^ oo 

by Hypothesis fiBll'). and the initial condition as well as the state constraints 
are satished by dehnition. Hence, by Lemma [21 there exist x{-) G a 

subsequence N" C N' and functions ^v(-) £ 1H^’^([0, r], R,'^), N G N", that 
coincide with yN{-) on [0, tArqr//ijvj] such that 

sup WyNit) — a;(t)|| —)■ 0 as N" 3 N ^ oo. 

tG[0,T] 

This contradicts assumption fflbll . and hence statement ca is correct. 

The second assertion is proved in a very similar way. Assume that r > r* 
and that there exist |/Ar(-) G ,/^Af, 0), N G N', for a subsequence 

N' C N. As above, we can extract another subsequence converging to a 
solution x{-) G -S'^(t), which contradicts = 0. □ 

5 Spatial discretization 

When the solution sets are computed in practice, it is necessary to discretize 
the underlying state space. Simple pathological examples show that in this 
situation, it is necessary to relax the state constraints in order to obtain 
convergence of the numerical scheme. The results obtained in this context 
as well as their proofs are similar to those in Section 01 

In addition to Hypotheses [^ and [HI we posit the following assumptions 
on the discretization of the state space. 

Hypothesis C. The spatial nets An G R'^ satisfy the following conditions: 

(1) For every x G R'^, there exists some y* E An with ||x — < Sn- 

(2) The relaxation parameter 6 n satisfies ^Ar/ZiAr—j-O as N ^ oo. 
Consider the discrete solution sets 

S'ffin, (3%, 6 n) := {yN e (3*, 6n) ■ yN,n e An, n = 0 ,..., i/} 
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with Q<u <N and 

+ (2 + Lhj\f)6]\f = L[P + l)/i^ + (2 + Lh]sf)Sisf. 

The following computations show that > 0 is a good choice for the blowup 
parameter. Hypothesis flC^) strengthens Hypothesis flB^) and ensures that 
(3* satishes Hypothesis flBl lh 

The next proposition implies that the fully discrete Euler trajectories are 
as close to the exact ones as the grid size allows. In particular, they do not 
cease to exist before the exact solutions violate the constraints. 

Proposition 8. For any x{-) G S'^{t,Xo) with r G (0,r*], there exists a 
sequence {y*)N with y*{-) G S'^^([r/hArJ,(57v) such that 

\\x{tN,n) - yti^N,n) |1 < <5iV, U = 0, . . . , [_T/h^l ■ (17) 

Proof. We construct the desired trajectory recursively. By Hypothesis flOl'). 
there exists an element 

G n H5 ^(x(0)) c n Bs^{Xo) n H 5 ^(A(o)) 

which satisfies the relaxed initial condition (E]) as well as the relaxed state con¬ 
straints dH]) and statement dn]). Assume that y^^,y% with {yt G 
S'3f{n,l3*,5N) and 

\\x{tNj-yt,k\\<^, k = 0,...n, 
have been constructed and that n < v. Then Lemma |2] yields 
dist(a;(fArn+i), + h^PitN^n, 2/v,n)) < (1 + LhN)hN + L{P + l)h^. (18) 
By Hypothesis dOl), there exists a point 

yN,n+l ^ n Bsj^{x{tN^n+l)) C A^r 0 Bsj^{A{tN,n+l)) 
satisfying the relaxed state constraint dS]). Because of ffTSj) . we have 
yN,n+l e yt,n + hNF{tN,n,yt,n) + ^/3#(0), 

SO that y*Q,..., y%n+i ^ trajectory of the fully discrete Euler scheme. □ 
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The following proposition shows that the opposite Hausdorff semi-distance 
between fully discrete and exact trajectories converges to zero. In particular, 
fully discrete Euler trajectories do not exist longer asymptotically than exact 
trajectories. 

Proposition 9. For any r G (0,r*], we have 


sup inf max \\yN{tN,n) - x{tN,n)\\ ^ 0 

as TV —)■ oo. For any r > r*, there exists E N such that 

S^f{[T/hN\,Pt,SN)=^ yN>Nr. 


Proof. The proof is the same as that of Proposition [71 provided the assump¬ 
tions of Lemma [3| can be verihed in the present setting. But for any sequence 
)n with G (3*,6isf), the linearly interpolated trajectories 

•) satisfy 


vti 


sup sup Il^^(t)|| < P + sup P^/hN < OO 

A''GNte(0,tjv, vgn 

as well as 

dist((f, 2 /^(f), graph(F)) 

< - {tN,it/h^\,yt{tN,it/hM\),ytitN,it/hM\))\\ 

+ dist ((f JV, p/hjv J) y* i^N, )) y* (tN, )), graph(F)) 

^ (-Zj -|- l)/i7v -|- 2.5]\fIh]\j' “h LSjsf —y 0 rs N —y oo^ 

and by the relaxed initial condition and state constraints, we have 

dist(|/^(0), Xo) < (jw —t 0 as N ^ OO, 
sup dist(|/^(fAr,„), A(tN,n)) < -t 0 as N ^ OO, 

n=0,...,N 

SO that Lemma [3] is applicable. □ 
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1=0.500000 


1=0.750000 




1=1.000000 1=1.500000 




1=1.750000 1=2.000000 




Figure 1: Approximate reachable sets for Example |5] with step-size h 


= 0.25. 
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t=0.500000 


t=1.700000 




t=2.500000 t=3.200000 




t=3.600000 t=4.000000 




Figure 2: Approximate reachable sets for Example |5] with step-size h 


= 0.05. 
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t=0.250000 


t=0.800000 




t=1.500000 t=2.200000 




t=2.800000 t=3.200000 




Figure 3: Approximate reachable sets for Example |5] with step-size h 


= 0 . 01 . 
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Hn 

Pn 

Sn 

T* 

0.2500 

0.3281 

0.0625 

OO 

0.2000 

0.2080 

0.0400 

OO 

0.1800 

0.1678 

0.0324 

20.1600 

0.1500 

0.1159 

0.0225 

10.3500 

0.1200 

0.0737 

0.0144 

7.0800 

0.1000 

0.0510 

0.0100 

6.0000 

0.0750 

0.0285 

0.0065 

4.8000 

0.0500 

0.0126 

0.0025 

4.1000 

0.0250 

0.0031 

0.0006 

3.5500 

0.0100 

0.0005 

0.0001 

3.2900 


Figure 4: Data for Example [5] with constants L = 1 and P = 2. Here, 
denotes the maximal interval of existence for the numerical trajectories. 




Figure 5: Visualized data for Example [5l 
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6 Numerical example 

We continue the simple, but instructive Example |5] from Section |H As it is 
difficult to visualize sets of trajectories, Figures [H [2] and [3] show the reachable 
sets 

= {yN,y ■ {yN,n)n =0 ^ S%* {v , , 5 n)} 

of the spatially discretized inflated Euler scheme for step-sizes h = 0.25, 
h = 0.05 and h = 0.01, i.e. the sets of all points that can be reached by 
trajectories of the numerical scheme at time t = 

In the hgures, the blue semiarc depicts the admissible set 

A = {x G 4- ^2 = 1, Xi < 0}, 

which is constant in time. The black line drawn around the semiarc represents 
the boundary of the admissible set A -|- Bsj^{0) of the numerical scheme. 
Blue dots are elements of the discrete reachable sets (3*, 6n), and the 

centers of the large red dots indicate the position of the exact solution at 
time t = Whenever dots exceed their corresponding admissible sets, 

this is due to limited plotting precision. 

The step-sizes have been chosen in such a way that we can observe the 
typical behavior of the inflated Euler scheme, as predicted by the analysis in 
the sections above. The exact solution is always approximated by a numerical 
trajectory up to the mesh-size of the spatial grid, which implies that some 
approximate trajectories live at least as long as the exact solutions. The 
approximate reachable sets, however, can be substantially larger than that 
of the exact dynamics. Moreover, numerical trajectories may live much longer 
than their exact counterparts. Figure [T] shows an extreme case in which the 
reachable set of the numerical scheme becomes stationary after time t = 1.5, 
and hence some trajectories are dehned on R_|_, while the exact solution 
ceases to exist at time t* = vr. In Figures [2] and [21 we see that the discrete 
reachable sets converge towards the true solution as hjv —)■ 0 and that the 
life span 

:= max{tAr,t. : S'}f{u,f3'^,5N) 7^ 0} 

of the numerical scheme converges to the exact value r* = tt as /iat —)■ 0 . 
This behavior is clearly recognizable in the data given in Figures 0] and |5l 
where the inflation and relaxation parameters (3^ and 6n a.s well as the life 
span are plotted as functions oi hj^. 
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7 Exploiting stability theorems 

Many papers in optimal control theory are concerned with the relationship 
between the solution sets 5'“(r) and for r G [0,T]. The specihcation 

of conditions, under which property fiDll) below holds, is subject of ongo¬ 
ing research. An overview over the relevant literature as well as instructive 
counterexamples can be found in [5]. 

It is certain that conditions, which ensure property fiDllj. are not satished 
by the class of examples we are aiming at, because property fiDll') implies 
S'^{t) ^ 0 for arbitrarily large r. Nevertheless, in this section, we will assume 
this property as given and describe briefly its consequences for the inflated 
Euler scheme. 

The ordinary Euler scheme (without inflation) has been analyzed in this 
context in |3], so that similar results should also be made available for the 
inflated scheme. We would also like to provide a template, how property (iDllj 
can be exploited to obtain short and simple proofs for linear convergence of 
a numerical scheme. 

To this end, we suppose the following hypothesis in addition to Hypothe¬ 
ses and [Bl 

Hypothesis D. The solution sets *S'“(r) and *S'^(t) enjoy a linear stability 
property, and the constraints are Lipschitz: 

(1) For any r G (0,r*], there exists some K = > 0 such that for any 

a:“(-) G 5'“(r), there exists some x^{-) G S^{t) with 

max \\x^(t) — x'^{t)\\ < K max dist(a;“(f), A(f)). 

te[0,r] te[0,r] 

(2) The state constraints are La-L ipschitz, i.e. there exists La> ^ with 

dist(A(fi), ^(^ 2 )) < LA\ti — ^ 2 ! Vfi, ^2 £ [0, T]. 

As Proposition [6] yields that the hrst Hausdorff semidistance between the 
exact and the numerical solution sets is zero, we only need to take care of 
the second semidistance. The proof of Proposition [TU] is remarkably short 
compared to the machinery set up in [3] to prove a similar estimate under 
similar assumptions. 
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Proposition 10. For any r G (0,r*], there exists C > 0 such that 

sup inf max \\yN,n - x(tN^ri)\\ < ChN- (19) 

yjvGSXr(L^/^JvJ,^iV,0) n=0,...,Lr//i;vJ 

Proof. A short computation shows that for arbitrary ypf G *S'^( [t//iatJ , (3n, 0), 
TV G N, the residual of the interpolated Euler trajectory yN{-) inserted into 
the differential inclusion ([T]) satishes 

dist(|/jv(t), F{t, yN{t))) < ‘2.L{P + Vt G [0, tw.Lr/hivj]- 

Therefore, Filippov’s theorem [H Theorem 2.4.1] guarantees that there exist 
Cl = Ci(L,P) > 0 independent of the choice of and a solution G 

5'“(tv,Lr/hivj) such that 


max |||/v(t) < Cih^. 

tG[0,tjv,[T/h^j] 

Since yN,n £ A{tN,n) holds for all n G {0,..., \j/h]^\}, Lemma [U and Hy¬ 
pothesis flD^j imply 

dist(|/jv(t),H(t)) < ||2/jv(t) - yN{tN,it/hj^\)\\ + dist{yN{tN^it/hj.,\), A{tN,it/hM\)) 

-hdist(H(tjv,Lt/^^j),H(t)) < {P + LA)hN, 

so that 

dist(a:)(r(t),H(t)) < C2hN Vt G [0,tN,[T/hj^\] 

with C 2 = Ci + P+La- By Hypothesis flDlH. there exist G S^{tN,[T/hN\) 

and a constant iV > 0 independent of the particular such that 

ll^iv(^) -a;v(^)ll < KC2hN Vt G 

and hence estimate flT^ follows with C = Ci + KC 2 - □ 

It is not difficult to see that the above result persists under a spatial 
discretization with parameter 6^ = hj.^. In this case, the hrst Hausdorff 
semidistance is still covered by Proposition |H1 
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8 Conclusion 


The aim of the present paper was to prove convergence of the inflated Euler 
scheme in the presence of state constraints and under weak assumptions, not 
excluding applications such as the collision avoidance problems discussed in 
m and ini. This aim is achieved, but the convergence statements are not as 
strong as one might wish, because we only showed pure convergence. It would 
be desirable to know if there were any conditions on the right-hand side and 
the state constraints that are substantially weaker than those imposed in [3], 
[5] and |6], but strong enough to specify a rate or a speed of convergence of 
some Euler-like scheme. In that sense, we leave a gap in the literature, which 
we are not able to fill at the moment. 
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